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Abstract 

A method for finding the world function of Robertson- Walker space- 
times is presented. It is applied to find the world function for the k = 
0,7 = 2, solution. The close point approximation for the Robertson- 
Walker world function is calculated upto fourth order. 

1 Introduction 

A fundamental invariant of a spacetime is the world function, which is given 
by 

<? = ±\t\ (l.l) 

where r is the geodesic distance, a is positive for timelike geodesies, zero for 
null geodesies, and negative for spacelike geodesies. This invariant contains 
a lot of information about an arbitrary Riemannian manifold. Once it is 
known the metric can be calculated from it using the coincidence limit of 
one point on the geodesic approaching the other, thus lim^.'^^ T a ^ — 9ab: 
and then the Riemann tensor can be calculated from the metric. Defining 
the Van Vleck determinant A = g 1 / 2 (x)det(a a b)g 1 ^ 2 (x') a Hadmard Green 
function is G(x,x') = z(87r 2 ) _1 A 1<2) (cr+ie) -1 , thus the world function, or an 
approximation to it, is needed when this Green function is required. These 
Green functions are used a lot in the theory of quantum fields in curved 
spaces, see for example Brown and Ottewill (1983) [||; and also they are 
used in the study of the motion of charged particles in curved spacetimes, 
Roberts (1989) §]. 

There are very few metrics for which the exact world function is known. 
It has been found for: DeSitter spacetime by Ruse (1930) Q|, for pp waves 
by Giinther (1965) Q, for Godel's solution by Buchdahl and Warner (1980) 
|||, and for the Bertotti-Robinson solution by John (1984,1989) @ @. Buch- 
dahl (1972) |1 has studies the production of the world function in Robertson- 
Walker spacetimes. In section | one new exact world function is found for 
the k = 0, 7 = 2 Robertson-walker solution, and previously known exact 
world functions for Robertson- Walker solutions collected together. For the 
majority of Robertson- Walker metrics it is impossible to find the world func- 
tion because intractable problems with elliptical integrals occur. For many 
applications it is only necessary to know the world function in the close 
point approximation and this is given in section [|. The conventions used 
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are: signature H ; a,b,c = 0, 1, 2, 3; k = 1, 2, 3; 

R%cd = W 6 - W + (i.2) 
= R a bda - 

The Robertson- Walker line element is 

ds 2 = N{t) 2 dt 2 - R(t) 2 {d X 2 + f(xf(d9 2 + sin 2 (6)dct> 2 )) (1.3) 
where R is the scale factor and 

!sin(x) if k = +1 

X if fc = (1.4) 

sinh(x) if k = —1 

in the conventions for fc of Hawking and Ellis (1973) ||p.l36. 

The lapse function N is arbitrary, and it can be absorbed into the line 
element; two common choices for N are N = 1, and N = R in which form 
the line element is conformal to the Einstein static universe. The Ricci 
tensor is given by 



R 



tt 



3R tt 3N jt R 



R NR 

?2 



g tj ( , RR tt 2Ri RR t N t \ . . 

Baa = zlL -2k f- f + ^ . (1.5 

3 R 2 \ N 2 N 2 N 3 J V ' 

Using Einstein's field equations the density \x and the pressure p for a perfect 
fluid are given by 



8vrG^ = — r + 



3A; 3R 2 i 



R 2 N 2 R 2 ' 

fc R 2 t 2N t R )t 2R, 



8irGp = — -77 + ^'^^ + (1.6) 

^ i? 2 N 2 R 2 N 3 R N 2 R V ; 

solutions to these equations with the equation of state 

p=( 7 -l)M, (1.7) 



have been given by Vajk (1969) [12|| . The speed of sound in a fluid is given 
by c dp/dfx, Hawking and Ellis (1973) [||p.91, thus the speed of sound for 
the equation of state L7 is 

c( 7 -l). (1.8) 
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2 The Geodesies of Robertson- Walker Spacetime 

The geodesic equations can be most elegantly derived by the Euler method, 
the Lagrangian is 

2£ = N 2 i 2 - R 2 g\fx i x j , (2.1) 

where: (i) 

(2.2) 

(ii) r is the geodesic distance, 

(hi) 2C = 1, 0, or — 1 for timelike, null, or spacelike geodescs respectively, 

(3) 

(iv) g\j is the metric of the three space of constant curvature, 
and i,j = 1,2,3. The momenta are 

P' = % = N 2 i (2.3) 

Pi KTJ — -tt Xi, 

OX 



which give 



dpt 
dr 

dpi 
dr 



= 2NN )t i 2 + N 2 i 

= 2RR t i±i + R 2 Xi . (2.4) 



The Lagrangian gives 



— = NN t i 2 - RRjfx^ 

|§ = -tfg&dtex'. (2.5) 



The Euler equation 

dC dp 1 

dx l dr ' 

after dividing by N 2 becomes 



(2.6) 



t+^ + ^9^ = 0. (2.7) 
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The Euler equation 



after dividing by R 2 becomes 



dC dpi 
dx l dr 



1R i 

X{ -\- — Hi 9kj,i'E X^ X — 0. 
K 



(2- 



(2.9) 



The geodesic distance u> on the three space of constant curvature is given 
by the well known expression 



cos{uj) = cos(x) cos (x') + sin(x)sin(x')cos(0) 



X 2 + x' 2 ~ 2xx'cos{@) 



if k = +1 

if k = (2.10) 



cosh{oj) = cosh(x)cosh(x') + sinh(x)sinh(x')cos(@) if k = — 1, 
where 



cos(Q) = cos(6)cos(6') + sin(6)sin(9')cos(4> - <j>'). 



Now 



,(3) 



g y'x l x 3 = I — 



doj \ dx % dxi 



Equations 2.1 and 2.12 give 



dr I doo du> 



2C = N 2 i 2 - R 2 Co 2 



and Equations 2.7 and 2.8 give 



Noting that 



Equation 2.9 becomes 



N f .9 RR f o 



(H}p j -r')., = (d;) ti = 0, 



ib + 2RR, t tto = 0. 



(2.11) 
(2.12) 

(2.13) 

(2.14) 

(2.15) 
(2.16) 



The geodesies are determined by 2.13, 2.14, and 2.16; 2.14 can be obtained 
by differentiating 2.13j and substituting 2.16, therefore only equations 2.13 



and |2,16| are used. Equation |2.16| can be integrated to give 



R 2 



(2.17) 
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where a is a constant independent of r. Rearranging 2.13 and substituting 



2.17 gives 



Ni = ±VR 2 Co 2 + 2C 



+ 2C 



(2.18) 



thus the dynamics of the geodesies are determined by |2.17 and 2.18. the 
integeral form of these equations are 



RN dt 



Va 2 + 21 gj R 2 ' 

aN dt 
RVa 2 + 2£#' 



(2.19) 
(2.20) 



takeing 2£ = +1 for timelike geodesies, these are just the equations de- 



scribing geodesies given by Robertson (1933) |1C]]. r is evaluated from |2.19 
with the positive sign and the limits of integration from T" to t and with a 



replaced by its value evaluated from 2.20. The the world function a is given 



from t by 1.1. 



3 Examples of exact World Functions for Robertson- 
Walker Spacetime. 



Example 1 The Einstein Static Universe. The metric is given by |l,3j with 

JV = 1, R = R , (3.1) 

where Ro is a constant. The ricci tensor is given by 

Rtt = 

Ikgij 



R; 



R 



(3.2) 



thus for / = the spacetime is flat. The pressure and the density are given 
by 

3k 

8ttGij = -24irGp = (3.3) 

Rn 



therefore the equation of state is [L7J with 7 = | . the world function is given 
by 

2a = (t - t') - R^uj 2 (3.4) 
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Example 2 The Generalized Milne Universe. For small time variation the 
metric can be expanded as a Taylor series in t 



N = 1 

R = R [l + H t - ^q H 2 t 2 = 0(t 3 )} (3.5) 

where Hq is the Hubble constant and go is the deceleration parameter. The 
qo ^ case proves to be intractable, so that only terms up to first order are 
considered 

N = 1 

R = R [l + H t}. (3.6) 

The Ricci tensor is 



Ru = 

(RlHj + k) 

thus when fc = — 1 and RqHq = 1 the spacetime is flat, and called the Milne 
Universe, which is just Minkowski spacetime with a point removed at the 
origin. The pressure and density are given by 

3(fc + R 2 H 2 ) 

8 ^ = - 2 ^ Gp =rWThW m 

therefore the equation of state is L5 with 7 = f • Gefining the new time 
coordinate 

t = -ET +t > (3 - 9) 
-no 

the metric is given by 

N = 1, 

R = R H t. (3.10) 
The world function is given by 

2a = t 2 + i' 2 - 2tt'cosh(R H uj). (3.11) 



7 



Transforming to the original coordinate system 3.£ this is 

2a = (t + t'f -4^J-+tj(l- + t'^j sinh 2 (^R H u;) (3.12) 

Example 3 DeSitter and Anti DeSitter Spacetimes. [] These are maxi- 
mally symmetric solutions to 

Rab = -Agab (3.13) 

with A > for DeSitter spacetime and A < for Anti DeSitter spacetime. 
For DeSitter spacetime the metric pressure and density are given by 

k = 0, N = R=^, X = \l'i, (3.14) 
8vrG/i = -8vrGp = A, 

and therefore the equation of state is given by |L7| with 7 = 0. For DeSitter 
spacetime the world funcion has been found by Ruse (1930) |11[| and is 



r = \f2~0 = — arc cos(Q), (3.15) 
A 



where 

g= 2^ ( * 2 + *' V) (3 ' 16) 
is the geometric distance in the five dimensional embedding space, and obeys 
the differential equation 

Q a Q a = A 2 (l-Q 2 ). (3.17) 

Transferring coordinates the result for DeSitter spacetime can be expressed 
as 



k = +l, N = l, R=\cosh{Xt), \ = J—, (3.18) 

A V 

cosh(r) = sinh(Xt)sinh(Xt') + cosh(Xt)cosh(Xt')cos(u) / X), 



footnote added 1999, a more complete treatment of the world function in DeSitter and 
Anti DeSitter spacetimes can be found in my work "Two Point Gravitational Energy" in 
preparation. 
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by symmetry the world function for Anti DeSitter spacetime can be found 
to be 



k = -l, N = l, R=\cos(Xt), \ = J^-, 

A y 3 

cos(r) = sin(\t)sin(\t') + cos(\t)cos(\t')cosh(u> / A). 



(3.19) 



Example 4- Robertson- Walker Spacetime with a Massless Scalar Field 
as Source. The world function for this spacetime is new and is derived here 
in detail, thus illustrating the method of Section |||. The metric is given by 



1.3 with 



N = R = bt 1/2 , 
where b is a constant. The Ricci tensor is given by 

-3 



(3.20) 



Rtt 
Rij 



The pressure and density are 

8vrG/i 

8nGp 



2t 2 ' 



bH 9ir 



_3_ / _L_ 

bH V + 4i 2 

— (h A 

bH V + 4t 2 



(3.21) 



(3.22) 



thus when k = o the equation of state is just 1.7 with 7 = 2, for which the 
speed of sound in the fluid equals the speed of light in a vacuum. 

Such a fluid is equivalent to a massless scalar field, see for example 
Roberts (1989) |J, with a Ricci tensor 



R 



ab 



for 3.21 <j) is given by 



-Vsin(t). 



The integrals |j| and |]20 become 
' bH dt 



2 r-„ --( 2a 2 

---Va 2 + bH [t — 

V Va 2 + bH 3 I b j 



li'j 



a dt 



2a 



Va 2 + bH b 2 



Va 2 + bH\l 



(3.23) 
(3.24) 

(3.25) 
(3.26) 
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Squaring |3.26 gives 



4a 2 



2a 2 - b 2 {t + if) = -2Va 2 + b 2 Wa 2 + b 2 t>, 



Squaring 3.27 and multiplying by orb gives the quadratic in a 



u A b A aWb 2 



16 



+ + « 4 ((t-t') -oj 2 ) = 



which has solution 



a 2 = \o~ x b 2 u?{t + t'± V&t' + uj 2 ), 



where 



2a e = (t - t') 2 - u 



then from 3.25 the world function is given by 



Y b ^=Yb T = ^e^it + t^^T^^ + tf' 2 



x(t- -a~ l uj 2 {t + t' ± VW +lo 2 )) 
-(l<r-W(t + t'± ViW + J 1 + 1') 1 ' 2 



x (t' - -aj V(i + t'± V&t' + UJ 2 )) 



(3.27) 

(3.28) 

(3.29) 
(3.30) 



(3.31) 



the metric is given by 1.3 with 



Example 5 The Tolman Universe. For the Tolman radiation Universe 

R = N = bt, (3.32) 
k = o, 



7= 3- 



The integrals |j| and ^20 can be evaluated to give 

at I bH a 2 , ,bt 

1 - — — arc sinh{ — 



a 2 2b 



a 



(3.33) 
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a . 6i, 
uo = — arc sinh{ — 
b a ' 



(3.34) 



and 3.34 cannot be inverted to give a. For the majority of solutions with 
equation of state L7 , the integrals 2.1S and 2.2C are elliptical and the prob- 
lem becomes intractable. The Tolman universe illustrates that the occur- 
rence of elliptical integrals is not the only problem that might arise; because 
in this case the integrals can be evaluated to give simple functions for r and 
u, however an expression for a cannot be extracted from these functions. 



4 The Close Point Approximation 

Introducing the notation 

ir ab - c = (x a - x' a ){x b - x' b ) . . . (x x c - x' c ), (4.1) 

and 

* n = <...<, (4.2) 
the world function can be expressed in the form 

2a = g ab 7r ab + h abc 7r abc + 0(7r i ). (4.3) 

Differentiating 

2a e = 2g ae ir a + Tr ab (g abye + h abe + h aeb + h eab ) + 0(tt 3 ), (4.4) 

thus 

°e<? e = g ab 7r ab + ±(g ab<c + 3h abc )n abc + 0{^), (4.5) 

now because 

2a = <r e a e i (4.6) 



Equations 4.3 and 4.5 can be equated to give 



1 

Kbc = -^9ab,c- (4.7) 

Carrying out the above approximation scheme to fourt order gives 

2<T = g ab TT - -g a b,cK 

+ j2 [^gab^-g'HlgaKe-gae^gcdj-gcf^^ + OiTr 5 )^) 
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For the line element 1.3 this gives 



2a = N 2 5t 2 - R 2 uj 2 - NNjSt 3 + RR t 5tu) 2 (4.9) 

1 /o*r2 , AATAT N, + 4 1 f D r> , RN ,t R A X+2,.2 



where 



-(3iVj, + 4NN tt )6t* - - \RR itt + — ^ j Sfu 

1 /^,A 2 w 4 + 0( ^5 )) 

St = t-t', (4.10) 



12 V JV 



and to is given by |2.11 
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